ESAIM: Probability and Statistics Will be set by the publisher 

URL; http://www.emath.fr/ps/ 



o 
o, 

(N| ; BEHAVIOR OF THE EULER SCHEME WITH DECREASING STEP IN A 

U ■ DEGENERATE SITUATION 



< 



c^ 



u 



Vincent Lemaire^ 



^^ ' Abstract. The aim of this paper is to study the behavior of the weighted empirical measures of the 

^ , decreasing step Euler scheme of a one-dimensional diffusion process having multiple invariant measures. 

• ' This situation can occur when the drift and the diffusion coefficient are vanish simultaneously. 

1 -^ I As a first step, we give a brief description of the Feller's classification of the one-dimensional process. 

j^ , We recall the concept of attractive and repulsive boundary point and introduce the concept of strongly 

repulsive point. That allows us to establish a classification of the ergodic behavior of the diffusion. We 

conclude this section by giving necessary and sufficient conditions on the nature of boundary points in 

terms of Lyapunov functions. 

In the second section we use this characterization to study the decreasing step Euler scheme. We 
give also an numerical example in higher dimension. 
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1. Introduction and framework 



Let / =]Z,'r[ denote an open (non-trivial) interval of the real line R. We consider the following stochastic 
differential equation 

dXt = b{Xt)dt + a{Xt)dBt, (1) 

S^ •■ where Xq is a random variable taking values in / and (^Bt).^„ a standard Brownian motion on R. We assume 

that b and a are continuous functions on / taking values in R, and that a is not degenerate on / i.e. Vx G /, 
(j'^ix) > 0. Then there exists a unique solution (Xt) , adapted to the completed Brownian filtration, such that 
1 1-^ Xt is continuous on [0, C[, where C — inf {t ^ 0, Xt — I or Xt = ?'} is the explosion time of the diffusion. 
In the first part, we establish a new "ergodic classification" for the process (Xt) in particular when 

C ~ +00. More precisely, we give the behavior of the sequence of empirical measures (t't),^„ — ( fn '^x ds ] 

according to the nature of the boundary points I and r. We characterize then the nature of the boundary points 
in terms of Lyapunov functions. The Lyapunov functions are usually used in high dimension, but there is a 
close link between these functions and the Feller's classification. This link makes it possible to more easily study 
the Euler scheme. 
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In the second part, we study the Euler scheme with decreasing step of a diffusion {X^^ on the real hne 
having (at least) a point A such that 5(A) = cr(A) = 0. In this situation we have 

Vx e] -oo,A[, P:r [^t e]-oo,A]] = I, and Va; e]A, +cx)[, P^; [Xt e [A, +cx)[] = I. 

In fact, the process {X^ has an ergodic behavior in I\ =] — cx), A[ or in I2 =]A, +oo[ according to the starting 
point X. But the Euler scheme is not continuous and may jump above the boundary point A. A legitimate 
question is then which are the weak limit of the empirical measures of the scheme ? We answer it in some cases. 

2. Results for the time continuous process 

We first introduce the scale function and the speed measure of the diffusion process {X^ solution of (^|. 
The next two sections are adapted from classical work on the Feller classification (see for instance [6], [5], [1] 
and [10]). 

2.1. Scale function and speed measure 

Definition 2.1 (Scale function). A scale function for the SDK J^l is defined for any c € / by 

VxG/, p{x)^ exp(- -iijdzjdj/. 

A scale function p is a strictly increasing function defined up to an afhne transformation. For the sake of 
simplicity, we call p the scale function of the process (X^)^^ 

We notice that the continuity of b and a, and the non-degeneracy of cr on / imply that p is in C^(/, R). 
Moreover, p is a solution of the following ordinary differential equation 

Vx e I, h{x)p'{x) + -a'^{x)p"{x) =0 i.e. Ap{x) = 0, 

and this property characterizes it. The probability that the process starting at x hits a point a Cz I before a 
point 6 S / is then expressed by using the scale function p. For any a £ I we denote Tq the hitting time of 
the one-point set {a} i.e. Ta = inf {i ^0, Xt — a}, and we consider a non-trivial interval ]a, 6[c / (strictly 
included). The function u{x) defined on ]a, b[ by Va; e]a. b[, u{x) — P^ [xTa^Tt — b] is solution to the system 

{Au = on ]a, b[ 
u{a) = and u{b) = 1, 

so that for all x £]a, b[, 

u{x) =P^[Tb<Ta\^ —- --. (2) 

p{b) ~p(a) 

This characterization is often used to define the scale function in a more general framework i.e. for continuous 
strongly Markovian processes which are regular in Dynkin's sense (Vx € /, Vy € /, Pa; [Ty < -l-oo] > 0) (cf. [11] 
or [10]). 

The following proposition gives another characterization of the scale function. 

Proposition 2.2. The process (p(x'l)] is a local martingale if and only if p is the scale function. 
Proof. For a proof in a more general framework, see proposition VII. 3. 5 in [10]. 
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If [p(xj) ] is a local martingale, then for all a < a; < 5 the process \p(xt "^ '') ) is a bounded martingale 
V /t:so V / t^o 

orem, we have 

p{x) = p{a)P, [Ta < Tb] +p(6)Px [Tb < TJ , 



and by the optional sampling theorem, we have 



which implies Q- 

If p is the scale function, then Ap — 0, and by the Ito's lemma applied to (a^t),>n with p we deduce that 

(p(xf)] is a local martingale. D 

The above proposition is very useful because it makes it possible to consider a one-dimensional diffusion as 
a Brownian local martingale up to a simple transform. 

Namely, the process (^t)j>Q defined for every t ^ by Yt — p{Xt) satisfies the following equation 

Yt^Yo+ [ .g(n)dB„ (3) 

Jo 

where Yq — p{Xo) e p{I) a.s. and g is defined by 

g(y) = /((^^')°^~')(y) ifyep(/), 

1 otherwise. 

The process (^t)j>Q can be seen as time-changed Brownian motion. The speed measure of (^t),>o evaluates 
how the change time affects the average time of exit from a bounded open interval. Let A be the generator 



of the process (^t)(>n and v the function defined on J —]a,b[c p{I) by v{y) = E^ 
inf {i ^0, Yt = z}. The function v is solution to the system 

{Av{y) = -1, a<y<b, 
\v{&) =:v{b) =0. 

Moreover by © we have Av{y) — \g'^{y)v"{y) for every y e]a, b[ and then 



TaATr 



where T, — 



(4) 



Vye]a,6[, v{y)^ I G^~^{y,z)^—-dz, (5) 

where G^ ^(y, z) is the Green function defined by G^ ^(y, z) — ^ '-, for every (y, z) G]a, 6p. 

Definition 2.3 (Speed measure). The speed measure of the time-changed Brownian motion {p{Xt))^^Q defined 
in is the measure M with density rfi — 2g^^ with respect to the Lebesgue measure. 

The speed measure of the process (Xt) is the image of M by p^^ and is a measure with density m = ^nrr 
with respect to the Lebesgue measure. 

As the speed measure of (^Xt)^^^ is the image of the speed measure of {Yt)^^^ by p^^ it follows from Q that 

r 1 /-P^'W 

\/ye]a,b[, v{y)^EyfaAf~^^ G, ^(y,p(z))m(z)dz. 

L J Jp-Ha) 
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Denoting a = p~^{d) and b = p^^{b), we have for every x s]a, b[ 

v{p{x)) ^ / Gp(^a),pib)ip{x),p{z))m{z)dz. 

J a 

Moreover p is a one-to-one function, then it is straightforward that £^(3.) Tp(^a) A Tp(^i,) — Ex [Ta A Tb]. Using 
^ we check that for aU ]a, b[c I and x e]a, b[ 

E, [Ta A n] = (1 - P. [n < Ta]) f (piy) - p{a))miy)dy + P, [U < Ta] f (p(6) - piy))miy)dy. (6) 

J a J X 

In addition the scale function and the speed measure provide a very useful expression for the infinitesimal 
generator A associated the SDE (Q. Indeed, we easily check that 

V/ e C2(/,R), Afix) = b{x)f{x) + la^x)nx) = -^ (^^) . (7) 

2.2. Feller classification 

The classification of one-dimensional diffusion process is due to Feller, in particular in [2] and [3]. In par- 
allel, the Russian school established similar results, but the two terminologies do not always coincide. For a 
comparison and a synthesis we refer to [6] . 

A first concept characterizing the behavior of the process (X^) solution of (^) in a neighborhood of a 
boundary point of/ =]/,r[ is the attractivity. 

Definition 2.4 (Attractivity). A boundary point A (A = / or A = r) is said to be attractive if lim \p{b)\ < +cx3. 

bei 

The function p is defined up to a strictly increasing affine transformation but the fact that the limit in A 
of p is finite (or infinite) does not depend on it. Similarly, since \p{x) — p{y)\ < -l-oo for all x,y and that p is 
strictly increasing we have 

I is attractive <=> Va; G /, lim(p(x) — p{b)) < +00, 
r is attractive <;=> Vx G /, lim(p(6) — p{x)) < -fcx). 

These equivalences are sometimes used to define the attractivity. We now show the following proposition which 
justifies the name of "attractive point" . 

Proposition 2.5. // A is an attractive boundary point, then for all a ^ I and all x in the open interval with 
endpoints a and A we have 

Px [Ta ^ Ta] > 0, 
with Ta defined by T\ = lim T^. 

6^A 

bei 

Proof. We detail the proof for the case a < a: < A. Firstly, since (ATt) is continuous we have by the mean- 
value theorem that the function (6 i— > T;,) (for 6 > a;) is strictly increasing. Thus Ta is the strictly increasing 
limit of Tb when b increases to A. Hence we have 



fl MTfa<T4-{rA<r4, 



be/nq 

a;<fc<A 
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and P^ [Ta sC Ta] = lim P^ [Tt < Ta]. By @ we deduce that 

Px Ta ^ Ta\ = lim — — -— . 

b^Ap{b)-p{a) 

However A is attractive, therefore hm (p(b) — p{a)] < +00, and the proof is complete. D 

Definition 2.6 (Repulsivity) . A boundary point A is said to be repulsive if it is not attractive, i.e. hm \p{b)\ — 

bei 

+CXD. 

Since p is strictly increasing and finite for any point of /, it is clear that 

I is repulsive <^ Va; G /, lini(p{x) — p{b)) ~ +00, 
r is repulsive <=> Vx E I, lim(p(6) — p{x)) = +cxd. 

b^r 

We also check by proposition 12.51 that A repulsive implies that for all a e / and x in the open interval of 
endpoints a and A, P^; [Ta > Ta] = 1. 

2.2.1. Attainability 

If A is an attractive boundary point, then a trajectory starting at a; G / hits A before another point a E I 
with strictly positive probability. But does this event occur in a finished time? Yes, if the point A is attainable. 

Definition 2.7 (Attainability). A boundary point A is said to be attainable if for all a G / and x in the open 
interval of endpoints a and A we have 

lim E^ [Tb A Ta] < +00. 

An attainable boundary point is attractive, and if A is attractive, then A is attainable if and only if 
Px [Ta < +00] > (cf. lemma 6.2 in [6]). 

These two concepts "attractivity" and "attainability" make it possible to determine the behavior of the 
diffusion in an neighborhood of a boundary point. Note that other concepts of the Feller's classification are not 
evoked here: regular point (reflective, absorbent, adhesive), exit point, natural point, entrance point. 

2.3. Behavior of empirical measures: the ergodic point of view 

Using Feller's classification, we can know the asymptotic behavior of one trajectory of the diffusion. But to 
establish the behavior of empirical measures and the recurrence of the process, the concept of attractivity is not 
precise enough. Indeed, several situations can occur: the boundary points +cx) and —00 are both repulsive for 
the Brownian motion and for the Ornstein-Uhlenbeck process. But for the Brownian motion it is null recurrent 
(in dimension one) and for the O.U. process it is positive recurrent. A new concept then is introduced: "strong 
repulsivity" . 

Definition 2.8. A repulsive boundary point A is said to be strongly repulsive if for every c € I we have 

\Jlm{y)dy\ < +00. 

A strongly repulsive point is a repulsive point such that the speed measure is finite in a neighborhood of A. 
We emphasize that this concept is defined from the repulsivity. Indeed an attractive boundary point A may 
satisfy \J. m{y)dy\ < +00 for every c € / (see the following example). 

Example 2.9. (1) Let / =]0, +oo[ and b and a continuous on /. Furthermore assume that for every 

x e [0, 1], b{x) = l^/x and a{x) = cx^l^ with c e]l, V2[. We have 



Va;e]0,l[, pl[x) =expf- / 



dy 1 = x" 



fP,yij1 
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and since c > 1, for every x g]0, 1[, p{x) = i_\t^2 X ~^ — i_\/^2 therefore lim p{x) = — -j-ir. The 



boundary point is thus attractive. 

Moreover, the speed measure is finite in a neighborhood of since 



1 ,.1 

m{y)dy = 



,2y3/2y-l/c^'^y ^2Ylc^-\ 



c-y/'^y- 



2 1 



(8) 



andce]l,y2[. 
(2) In the case of the Ornstein-Uhlenbeck process defined on R by 



dXt 



1 



'-XtAt + dBt, Xo^xgR, 



(9) 



the speed measure mix) Ax is the Gaussian probabihty, and the boundary points — oo and +oo are thus 
strongly repulsive. 

We recall now the main ergodic result for one dimensional diffusion process {Xt)^ . A process (Xt)^ is 
said to be recurrent if for all a and b in / we have Pq \Tb < +oo] = 1. Moreover the process is called positive 
recurrent if Ea [Tb] < +oo and Efc [Ta] < +oo, and null recurrent otherwise. 

Theorem 2.10. We suppose that (Xt) (solution of ^ and with speed measure m) is recurrent on I. If f 
and g are two non negative measurable functions such that 



{f{x) + g{x)) m{dx) < +oo, / g{x)m{dx) ^ 0. 
/ Ji 



Then 



f{Xs)ds / f{x)m{dx) 



liX,)ds 



g{x)m{dx) 



Proof. For a detailed proof we refer the reader to [11] or [4]. 

In the sequel we will denote by (^'t)^^(, ^^^ empirical measures of the diffusion, i.e. 



D 



1 /•* 
yt > 0, vtidx) = - / ,5x,(dx)ds. 
t Jo 



By the above theorem and the concept of strongly attractive boundary point we establish the following 
classification of the ergodic behavior of the diffusion. 

Theorem 2.11. We recall that C = inf {t ^0,Xt^l or Xt = r] . Then 

• if I is attractive and r is repulsive then xl ——> I, 

• if I and r are attractive then 



lim xj — I 

t — i-+oo 



IP 



lim XI = 

t — > + C)0 



pir-)-piXo) 
p{r-) -p{l+) ■ 



if I and r are repulsive then the diffusion is recurrent and does not explode (Q — +oo a.s.). More 
precisely 
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- if I and r are strongly repulsive (i.e. the speed measure is finite) then the diffusion is positive 
recurrent and vt ^ v a.s. where v is the normalized speed measure. Moreover 

yfeL\,y), I f f{X,)ds^ [ fixMdx) a.s. 
t Ja Jr 

- if I is strongly repulsive and r is (simply) repulsive then the diffusion is null recurrent and if the 
empirical measures are tight we have 

1 /•* 

- / (5xsds ^ Sr a.s. 

- if I and r are not strongly repulsive then the diffusion is null recurrent, and if the empirical 
measures are tight then any weak limit 0/(7/0 '^-^s'^* ) '■^ "^ measure with support {/, r}. 

We first prove the following lemma. 

Lemma 2.12. // the two boundary points I and r are repulsive then the diffusion is positive recurrent if and 
only if its speed measure is finite. 

Proof. By definition, the diffusion is positive recurrent if and only if for all a and b in /, Ea [T^] < +00 and 
E{, [Ta] < +00. Let I < a < b < r. By symmetry it is sufficient to prove that 



E„ [Tb] < +00 ^ / 



m{y)dy < +00. (10) 



Firstly, I is repulsive and T; = lim^^i T^ therefore Ea [!&] = Ea [U A Ti] = hm^^, Ea [!& A T^]. Moreover by 
(O we have Vx £]l, a[, 

Ea [n A rj = Pa [n < T,] / {p{b) - p{y)) m{y)dy + Pa [T, ^ n] / {p{y) - p{x)) m{y)dy, 

J a J X 

and since / (j>if>) — Piu)) ''TT-{y)dy is finite and does not depend on x, the limit when x tends to I of Ea [Ti, A T^] 
is finite if and only if 

lim f Pa [T^ s; Tb] / {p{y) - p{x))m{y)dy j < +00. 
As Pa [T. ^ Tb] = ggEfg we have 

Pa [T. < Tb] r {p{y) - p{x))m{y)dy = {p{b) - p{a)) f ^M^l^m(y)dy, 

Jx Jx P(0)-P[X) 

= ipih)-p{a)) I Vy [Tb < T,] m(y)dy, 

and it follows that Ea [Tb] < +00 if and only if lim / Py [Tb < T^] m{y)dy < +00. Since T^ strictly increases 

to T;, Pj, [Tb < Tx] increases to Py [Tb < T/] = 1 because / is repulsive. The monotone convergence theorem 
yields ^. D 

Proof of Theorem \2.11\ The first two items are proved in [5] (Proposition 5.22). We recall the proof of the first 
item. 
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• Suppose that I is attractive and r is repulsive. By definition of the scale function we have for all 
l<a<x<b<r 



inf Xt ^ a 

OsSt<C 



^ Px [xTaAn = «] 



p{b) -p{x) 
p{b) ~p{a)' 



Increasing 6 to r we obtain P^; [info^t<^ Xt ^ a] — 1 since r is repulsive. The limit when a decreases to 
I also gives 



On the other hand 



sup Xt — b 



inf Xt = I 



^ Px [xTi/\n 



(11) 



b\ = hm —— -— . 

a^i p[b) — p{a) 



and taking the limit when b increases to r we obtain supg^j^^ Xt < r a.s. It remains to prove that the 
process (a^^) . is almost surely convergent. Since [p{x{) ) is a local martingale and I is an attractive 

boundary point, the process {p(x:) ~\mia^ip{a.)) is a positive continuous local martingale. By 

Fatou's lemma it is a positive continuous super-martingale which is then almost surely convergent. We 

conclude using the continuity oi p~^. 

If I and r are repulsive then in the same way that we obtain Hll() we have 

inf Xt — I a.s. and sup Xt = r a.s. 

The diffusion is thus recurrent on / and C = +cxd a.s. Moreover, by Lemma 12.121 we know that the 
recurrence is positive if and only if the speed measure is finite. 

- If the two boundary points are repulsive then the speed measure is finite and by Theorem l2.1()l 
we have 



V/ 



GL^M, i / f{Xs)ds^ f fdv a.. 

^ Jo JR 



where v is the normalized speed measure. 

- If / is strongly repulsive and r is repulsive then the diffusion is null recurrent. Considering an 
increasing sequence of continuous functions with compact support (<?n(a;))„>]^ such that g„(x) — > 1 and 
Vn ^ 1, ^^gn{x)ra{Ax) ^ 0, we obtain by Theorem 12. 101 



V/eLi(m), - f f{Xs)ds 
^ Jo 



0. 



(12) 



On the other hand, we consider a sub-sequence (^'a(t))f>n '^^ i^t)t>ri converging to a measure v (the 
empirical measures are tight). Let / be a continuous function with compact support such that supp(/) C 
[/, r[. As i^a{t) =^ J^ we have 

f{Xs)ds ^ / /d^. 



a(t) Jo 

The boundary point I is strongly repulsive and supp(/) C [I, r[, thus / is integrable with respect to m 
Hence (|12|) implies / fdiy ~ 0. The interval [/, r[ satisfies 



V/ e Ce(/), supp(/) C [Irl => v{f) = 0, 
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therefore supp(t^) = {r}. Since v is normalized we have v = 5r- 

- In the same way, if the two boundary points are strongly repulsive then any weak limit of \t^t) f^n. 
is a measure with support {l,r}. 

D 

2.4. Attractivity and Lyapunov function 

In order to study the behavior of the Euler scheme (with decreasing step) we establish a link between the 
concepts of attractivity, repulsivity and strong repulsivity, and the Lyapunov functions. Let us notice that there 
is few work which relates to this subject. Indeed, the Lyapunov functions are useful in high dimension and the 
Feller's classification is established for one-dimensional processes. 

In the sequel, we will denote by Ja C / an open (non-trivial) interval included in / with endpoint A. 

Proposition 2.13. Let A be a boundary point (finite or infinite, lefi, endpoint or right endpoint) of I . The 
following statements are equivalents 

(1) A is a repulsive boundary point of I if and only if there exists a neighborhood Ja d I of A and a strictly 
monotone function V G C^(Ja,R'+) such that 

lim V{x) = +00 and \/x E Ja, AV{x) s^ 0. 

(2) A is a strongly repulsive boundary point of I if and only if there exists a neighborhood Ja <Z I of A and 
a strictly monotone function V G C^(Ja,R+) such that 

3e > 0, Vx e Ja, AV{x) «: -e. 

(3) A is a attractive point of I if and only if there exists a neighborhood Ja d I of A and a strictly monotone 
function V G C^( Ja, R+) such that 

sup V{x) = V{A) < 4-00 and Vx e Ja, AV{x) ^ 0. 

xGJa 

Proof. We give the proof when A is a right endpoint of /. Then Ja is an interval ]c, A[ with c G /. 

• - We suppose that there exists a neighborhood Ja of A and a function V G C^(Ja,R+) such that 
limj,^A y{x) = +00 and AV ^ on Ja. For every x E Ja we have 

m[x) \ p (x) J 

hence V /p' is decreasing on Ja. There also exists C > such that V'{x) ^ Cp'{x) for every x g]c, A[ 
since p' > 0. It follows that lim p{x) — p{c) — +00 because V tends to infinity in A. 

x~^A 

- Conversely we must find the good Lyapunov function V. Let c > be such that p{c) > (c exists 
because A is repulsive). Since p is strictly increasing we have for every x g]c, A[, p{x) > p{c) > 0. We 
also define the function y on ]c, A[ by 

Va;G]c, A[, V (x) ^ p{x) - p{c) . 

The point A is repulsive thus V increases to infinity when x tends to A. Moreover V G C^Qc, A[, R+) 
and AV = 0. 

• - Let V G C^(Ja, R+) be such that lim^^A V{x) — -l-oo and £ > such that AV ^ — e. Thus we have 



AV {x)m{x)dx ^ —e / m{x)dx, (13) 

Ja -^ J A 
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1 fV'{x)\' 
and since AV{x) — — —^ — -— ) we obtain for every c G Ja, 
m[x) \ p (x) ) 



^ V'{x\ V'{c\ 

AV{x)m{x)d.x = lim -^ - -^. (14) 



s-tA p'(a;) p'(c) 
By p3f) we derive that 

1 (y\c) ,. V'{x) 



I Tn(a:)dx ^ - ,, , lim 

Jc £ KP'ic) --^ 



e \p'{c) x^A p'{x) J 
As the functions V and p are increasing on Ja we have hm — — — - ^ 0, which gives f m(x)dx ^ C 

x^A p'[x) 

{i.e. A strongly repulsive). 

- Conversely we assume that A is strongly repulsive. Let c G / and V the function defined on ]c, A[ by 

Vxe]c, A[, V{x) ^ \p'{y) m{z)dz\dy. 

It is clear that V G C^Qc, A[, R+) and that for every x g]c, A[, V'{x) — p'{x) J m{z)dz. Moreover 

V.Te]2;,A[, AV{x) = -l 

• We prove (3) in the same manner as (1). For the converse we consider the function V{x) = p{x) — p{c) 
on ]c, A[ with c such that p{c) > 0. 

D 

The above Proposition provides a useful criterion to know the nature of a boundary point. However, when 
the boundary point A is finite, a "natural" Lyapunov function has a minimum at A and it is not the case of 
V. But an easy transform allows us to obtain this property. This is the interest of the following Corollary. 

Corollary 2.14. Let A a boundary point of I. 

(1) A is a repulsive boundary point of I if and only if there exists a neighborhood J a d I of A and a strictly 
monotone function v € C^(Ja,R+) satisfying v{A) — 0, such that 

VxSJa, Av{x)^laHx)^-^^^^, (15) 

2 v(x) 

(2) A is a strongly repulsive boundary point of I if and only if there exists a neighborhood J a <Z I of A and 
a strictly monotone function v £ C'^(Ja,R+) having a minimum at A, such that 

3e > 0, Vx e Ja, Av(x) ^ ^^^(rE)^^^^^ + ^^(x), (16) 

2 v{x) 

(3) A is an attractive boundary point of I if and only if there exists a neighborhood J a C / o/ A and a 
strictly monotone function v G C^(Ja,R+) having a minimum at A, such that 

VxeJa, Av{x) < -a^xy-^^^P^. (17) 

2 v(x) 
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Proof. Let Ja a neighborhood of A strictly included in /. We consider the case in which A is the right endpoint 
of / i.e. Ja =]c, A[ with c G I. We define first for L ^ the function (J>l by 

0L:[O,exp(L)[^R+, 

X t-^ — hi{x) + L. 

The function (p^ is a strictly decreasing one-to-one C°° function. Moreover, for every C^ function v with values 
in [0,exp(L)[ we have 

^(^,o.) = -:^ + Vi^. (18) 

V Z v^ 

• - We assume that there exists v e C^([c, A],R+) strictly monotone satisfying u(A) = and (|15|) . We 
also define on ]c, A[ the function V by Va; €]c, A[, V{x) = {(pL o v) {x) with L = ln(7;(c)). It is a strictly 
monotone function which tends to infinity in A. By (|15() and (|18|l we obtain AV ^ on ]c, A[. The 
proposition (|2.13|) implies that A is repulsive. 

- Conversely if A is repulsive then there exists V £ C^Qc, A[, R+) strictly monotone which goes to +cx3 
when X tends to A. We also define v — cj)'^ o y on ]c, A[ with L = ij^ixe]c,A[V{^)i and we extend if 

by continuity on [c, A] letting v{c) = 1 and v{A) = 0. By AV ^ and H18I) we have Av ^ ^a'^-^^^-j— on 
]c,A[. 

• - Let V a strictly monotone function on [c, A] such that v{A) = 0. For 9^0 and L = ln(ti(c) + 9) 
we define V{x) = {4'L,e ° 'v){x) for every x g]c, A[. This function V is strictly monotone on ]c, A[ and 
admits a limit (finite or not) when x tends to A. From (|16|) and H18|) we deduce that 

Vxe]c,A[, AV{x)ii-e. 

- Conversely we consider the function v ~ (p^ o V^ on ]c, A[ with L — inixe]c,A[ ^i^) and we extend if 
by continuity letting v{c) = 1 and v{A) — lima;^A exp(— T^(a;)). 

• The proof is similar to the two firsts items. 

D 

Remark 2.15. If A is such that 6(A) — cr(A) = 0, then A is a critical point for the equation u' = h{u). It is 
worth noting that there exists a link between the nature of the critical point A for the ODE u' = b{u) and the 
nature of the boundary point A for the SDE. 

Indeed, if A is a stable critical point then there exists a Lyapunov function F E C^ such that F'b{u) < 
for every u in a neighborhood of A. If F' /F is decreasing, the above corollary implies that A is an attractive 
boundary point for the SDE. 

If A is an unstable point for the ODE u' — 6(u), it may be repulsive, strongly repulsive or attractive for the 
SDE, as shown in the following example. 

Example 2.16. We consider, like in [9], the function y : R ^ R+ defined by 

nx) = |^f"^'^"?^^' '^'"'^^' and 6(x) = |^'/V^r(")) ifN^3, 

and b = —V. The ordinary differential equation u' = b(u) has 3 critical points: —3, and 3. The points —3 
and 3 are stable and the point is unstable. 

Let c €]0, 2[ a parameter and a defined by (t{x) — ex. We consider the process (^Xt).^„ solution of the SDE 

dXt — b{Xt)dt + a{Xt)dBt. It is clear that the point is a boundary point for (Xt) . Moreover we check 
that 

AV{x) = - (4 ^ c^) x^ + 12 sgn(a;) if |a;| ^ 3. 



12 
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Figure 1. Fonction V 



Figure 2. Drift b ^ -V 



By Proposition 12 . 1 jl the points — oo and +cx) are then strongly repulsive. 

On the other hand, we use Corollary 12.141 to determine the nature of the boundary point according to c. 
Assume that / =]0, -|-oo[ and let v the function defined on [0, -|-oo[ by v{x) = x^. We have 



Vxe]0,3[, Av{x) ^ {I + c^)x^ 



9 



and -a'^{x)- — t-t— 
2 ^ ' v{x) 



= 2c^x'^, 



(19) 



and then 

• if c < 1 the boundary point is strongly repulsive (for ]0, -|-cx3[ and ] — oo. 0[ by symmetry). Indeed the 



condition ((TB|l is satisfied with e = ^-^ and Ja = ]0, 3i/ ^^^ [. 

• if c > 1 it is easy to check that the boundary point is attractive. 

• if c = 1 we consider the function v{x) — a;exp(x) and we check that is a repulsive boundary point. 

Thus the nature of the boundary point (which is always stable for the ODE u' = b{u)) may change according 
to c. By Theorem 12 . 1 1 1 the ergodic behavior of {XtY^ is the following: 



• if ce]l,2[ thcnXt 

rt 



■ (for every starting point Xq), 

if c = 1 then j J^ Sx^ {dx) ^ ^Oj 
if c < 1 then for every / S L}{m) 



' J fdv- ifXoe]0,+oo[ 
fiX,)ds ^^ { /(O) if Xo - 

Jfdiy+ ifXoe]-oo,0[ 



with i^_ the invariant measure of {Xt)^y„ on ] — oo, 0[ and j/-|_ the invariant measure on ]0, -foo[. 



3. Behavior of the Euler scheme with decreasing step 

We consider now the Euler scheme (Xn) >q built using a positive sequence (7„) ^ going to 0. We assume 
that (7„) ^ satisfies lim„ X]fe=i 7fc — +o^ ^^d we denote r„ = X]fe=i Ik- The inhomogeneous Markov chain 
(^")„^o ^^ defined by 

Xn+l — Xn + "fn+lK^n) + \/jn+lCr{Xn)Un+l, 
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with (Un) >, a real white noise i.e. a sequence of i.i.d. random variables such that E [Ui] ~ and var([/i) = 1. 
Furthermore, we assume that Ui is a generalized Gaussian (cf. [12]) i.e. such that 

3k > 0, V6' e R, E[cxp(6lL/i)] s$ exp j ^^ j . 

For example Ui is a standard Gaussian or a Bernoulli random variable. A consequence of the generalized 
Gaussian property is the following 

Va^O, P[\Ui\^a]^cxp(-^Y (20) 

We consider a diffusion (Xt)^^„ on the real line with b and a continuous on R. Moreover we assume that b 
and a have sublinear growth i.e. 

3Cb>0, \bf ^Cb{l + \x\), and 3^ > 0, \af ^ 0^1 + \x\). (21) 

In the sequel A denotes a finite point of R such that 6(A) = o'(A) = 0, and Ja denotes an open interval with 
endpoint A (Ja =]A, A + e[ if A is a left endpoint of / =]A, +oo[ and Ja =]A — e, A[ if it is a right endpoint 
of /=] -oo,A[). 

3.1. Euler scheme 

We prove the following Theorem which gives the behavior of one trajectory of the Euler scheme \Xn) >„ in 
this degenerate situation. 

Theorem 3.1. We assume that a satisfies (t{x) ^ for every x g] — oo, A[U]A, +oo[ and that there exists 
Ua =] a — e, a + e[ with e > 0, and a convex function v G C^{Ua, R.+) satisfying v(A) = and 

yxeUA, {v'b){x)^0 and 3c^ > 0, Vx e Z^A, \{v'a){x)\ r^ c^v{x). (22) 

// the step sequence (7n)„>]^ satisfies VC > 0, X]n>i ^^P ( ~ ) ^ "'"°° then the Euler scheme jump above A a 
finite number of times i.e. 

P[3no ^0, Vn^no, X„e]-oo,A[] +P[3no ^0, Vn ^ no, X„g]A,+oo[] = 1. 

We first prove the following lemma. 
Lemma 3.2. We assume there exists a convex function v G C^(JajR'+) satisfying v(A) — and 

Vx e Ja, {v'b)ix) ^ and 3c„ > 0, Vx e Ja, \{v'a){x)\ ^ c„v{x). (23) 

Then, on the event {Xn G Ja} 

P[A e {Xn,Xr,+ l)\Tn] Sj CXp {-^ 

Proof. We assume that A is a left endpoint of / and we denote by A„+i the event {A e (X„,X„_|_i)} (the 
geometric segment with endpoints Xn and X„_|_i). Since v is continuous and f (A) = it is clear that 

An+l = {3t e [0, 1], Xn + t{Xn+l - Xn)) = A} , 

C {3t e [0, 1], v{Xn + t{Xn+l - Xn)) = 0} . (24) 
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Consider t E [0, 1] such that A = X„ + t{Xn+i — Xn). As v is C^ on Ja, Taylor's formula gives 



viXn + t{X„+i - X„)) = viX„) + v'{Xn)t{X^+i - X„) + 



«"(en+l) 



t \Xn+l — Xn) , 



with ^n+1 e]A, X„[- The convexity of w implies 

= v{Xn + i(X„+i - Xn)) ^ v{Xn) + t-1n+l{v'b){Xn) + t,J^i^{v' a){Xn)Un+l. 

Since v'b ^ on Ja we have from H24|l 

An+i n {X„ G Ja} c {3f e [0, 1], v{Xn) + tV7^(v'^)(^«)[/n+i s; 0} n {X„ e Ja} 
c {3i e [0, 1], iV7;;Tr|(fV)(x„)c/„+i| ^ i-(x„)} n {x„ e Ja} . 



Hence for any n ^ 0, we have on the event {Xn G Ja} 

PK+i \Tn] S=;P |C/„+l| ^ 

s;p |c/„+i| ^ 



i;(X„ 



,/i;;ri\W<T){Xn)\ 
1 



-' rz 



'^ n 



by the domination assumption (|23|) on v'a. We conclude using property (|20|l of the random variable Ui 
Proof. By Lemma 13.21 we prove easily that for every n ^ 



P[A e {Xn,Xn+l)\Tn] ^ exp ( - ^^— 



on Xn e Wa- 



D 



(25) 



We now consider the event {Xn ^ Z-/a}- Then we have 



{A e (X„,X„+i)} = 3i e [0, 1], Un+i = 



A-Xn 



t^/ln+l(^{Xn 



C <^ \Un+l\ ^ 



\A-Xn 



V7n+l|cr(X„)| 



V7n+1 



- V7«+i 

\biXn)\ 
W{Xn)\ 



bjXn) 
<T{Xn) 



As the drift b is dominated by Cb(l + \x\), we have 



{Ae{Xn,Xn+l)}c\\Un+l\^ 



|A-X„| 



/7;:^Cb(i + |x„|) 



V7n+1 



Cbil + \Xn\) 
HXn)\ 



and using the triangular inequality and |X„ — A| ^ e we prove that '. .y "J J^ — iT^ 
there exists ni ^ and C > such that for any n ^ ni, 

P [{A e (x„, x„+i)} n {x„ ^ Wa} I Tn] < P I 

<P ||f7«+i| ^ 

using Icrl < yC7^%/y. 



We also deduce that 



\Un+i\ > _^ ^f + 1-^-1) ^ n {X„ ^ Wa} 



V7n+i k(A:„)| 



.F„ 



Cay^, 



n+1 



n {x„ ^ Ua} 



Tn 



(26) 
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From l|^ and ll^ combined with it^ . we get 

C 



3ni > 0, 3C > 0, Vn ^ m, P[A e (X„,X„+i)|j^„] s^ exp 



7«+i 



By the condition on the sequence (7„) . and the conditional Borcl-Cantelh Lemma we deduce that the event 
{A G {XmXn+i)} occurs a finite number of times. D 

Remark 3.3. The condition on the step sequence (7,1) ^ is not restrictive. Indeed, it is satisfied for (7„) ^ 
defined by 7„ — jon~^ with 70 > and r s]0, 1], or 7„ = log(n)~'' with r > 1. 

The technical assumption "?; convex" is not very restrictive in practice. The important point to note is the 
condition v'b ^ which implies that A is unstable for the ODE u' = b{u). But A may be repulsive as well as 
attractive for the SDE (cf. Remark |2. 15(1 . The condition v'cr — 0{v) in a neighborhood of A is very important 
and it seems difficult to relax it. 

Remark 3.4. It is easy to extend the above Theorem to the case of finitely many boundary points A^. If 
for every Aj, there exists a neighborhood Wa; and a convex function Vi G C^(Z//Ai, R-) such that v(A) = and 
satisfying (|22ll . then 

P[3ie{0,...,l}, 3no>0, Vn^no, X„g]A,,A,+iO = 1, 
with Aq = —00 and A;+i = +00. 

3.2. Weighted empirical measures 

Let (r]n) >-, a positive sequence, called weight sequence, such that iJ„ = Sfc=i Vk increases to +00 when n 
tends to +00. We define the weighted empirical measures (i^^)„>;^ by 



1 " 



'^ fe=i 
In this section we assume that the diffusion satisfies a stability condition i.e. 

3a > 0, V|a;| > Af, xb{x) + -a'^ {x) i^ - ax^ . 

This condition implies that the points — cxd and +00 are strongly repulsive and that the empirical measures 
iyt) f^r. are tight. Since b and a have sublinear growth, this condition implies also the tightness of the weighted 
empirical measures (i^^) . of the scheme and that any weak limit is an invariant probability for the diffusion 
(cf. [7] or [8]). 

A consequence of Theorem 13.11 is the following proposition which describe the convergence of {y^) y, ac- 
cording to the behavior of b and a in a neighborhood of A. For more clearness, we parametrize b and a. 

Proposition 3.5. Let A the unique point ofR such that &(A) — a{A) = 0. We assume that in a neighborhood 
of A we have b{x) — sgn(x — A)pi,{x) and a{x) = Pa{x) with pi, ^ 0, 

Pb{x) ^ Cb\x - Af and a{x) ^ c„\x - A^ (27) 

where (3, <i, Cb and c^ are positive real numbers and <; ^ 1. // the step sequence (7„) ^. satisfies VC > 0, 
J2n^i exp(-C/7„) < +00, then 

• if 1 + P — 2i; > then A is an attractive boundary point and J^^ => (5a , 

• if 1 + P — 2i, = and Ca > V2cb then A is an attractive boundary point and i^^ ^ ^a, 
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• if 1 + (3 ~ 2<^ — 0, Ca < ^/2cb and (3 ~ I (which implies <; ^ 1) then A is a strongly repulsive boundary 
point and v^^ =^ v^ or i>''^ => V- , 

• z/ 1 + /3 — 2<j < 0, Co- ^ \/2ct, and (3 e]0, 1] then A is a strongly repulsive boundary point and v^ => v^ 
orvl ^ v_, 

where Vj^ is the invariant probability on] A, +oo[ and v^ the invariant probability on ] — oo, A[. 

Proof. To simplify notation, we assume without loss of generality that A = 0. 

Let Z^A — ]^ — e, A + e[ a neighborhood of A and the convex function v{x) = x^ . Then for every x G Z^a we 
have 

v'h = 2x Bgx\{x) pi,{x) ^ 0. 

Moreover v'a ^ Ca\x\'' with <^ ^ 1 hence there exists C > such that \v'a\ ^ Cv. By Theorem 13. II we know 
then that the scheme lives in ] — cx3, A[ or in ]A, +cx)[ after an almost-surely finite random time. 
Furthermore, with vix) = x^ we have 

Vx G Z^A, Av(x) = 2\x\pb{x) + u'ix) = 2c6|a;|'+^ + c^ |x|'' + o(|xp+'^)''('^)) , (28) 

and \a''{x)^^^^^ = 2a\x) ^ 2cl\x\^\ 
2 v[x) 

- If 1 + /3 > 2^, there exists a neighborhood of in which Av < ^cr-^^^^. Hence by Corollary 12.141 is 
an attractive boundary point. The sequence of weighted empirical measures (f^j . of the scheme is tight on 
[0, +cxd[ and on ] — oo,0], and any weak limit is an invariant probability. However, So (the Dirac at 0) is the 
unique invariant probability on [0, +oo[ or on ] — cx3, 0]. Thus any weak limit of (i^n) >i i^ ^Oi which proves the 
first item. 

- If 1 + /3 = 2<;, Av{x) = {2cb + cj)|a;| *■ + o( |a;| *" ) . If Cg- > \/2cb there exists a neighborhood of in which 

Av < 2ct^ and we conclude as above. 

If /3 = 1 and c^ < \/2cb, we have for every x in Ua, Av{x) ^ 2(t^(.t) + ex'^. By Corollarv l2.14l the point 
is then strongly repulsive. Any weak limit of (z^^) . is a probability on ] — cxi,0[ or on ]0,+oo[, therefore we 
have t'^ => i/_ or v^ => v^ (we recall that the boundary points — oo and +oo are strongly repulsive). 

- The proof for the case 1 + /3 < 2<^, c^ ^ \/2c6 and /3 e]0, 1] is similar. D 

To illustrate this result we take the same example as in the previous section (Example 12. lt)|) . 
Example 3.6. Let b and a be defined by 

, , , I —2(x — 3 sgn(a;)) if |a;| ^ 3, , , , . , ,„ „r 

b(x) = I } ^ & \ JJ I I ^ ' and a(x) = ex with c e]0,2 . 

\-j^.T-^ + i:r if |x| «:3 ' ^ ' ^ ^ 

We recall that the nature of the boundary point for the diffusion {Xt)^^^ depends on parameter c. If c g]1, 2[, 
is an attractive boundary point, if c = 1 it is repulsive and if c £]0, 1[ it is strongly repulsive. 

By the above Proposition, we know that the weighted empirical measure (I'n) >i weakly converge to Sq when 
c g]1, 2[ and to i^+ or :/_ when c g]0, 1[. Note that the convergence to v^ or ;/_ does not depend on the initial 
condition and is not previsible. 

We give a representation of the density of v approximated by f^ with n = 10^. More precisely, we have 
discretized the interval [—2,8] using 200 intervals li of length 0.05 and we have computed t'^(l/J for each li 
with n = 10^. The step sequence (7^) ->-, is defined by 7„ = n~^'^ and the weight sequence yrjnj ^, is defined 
by rjn ~ I. The results of this approximation of the stationary density are given in Figure El for different values 
of c. 

We remark that for a small noise (c — 0.1), the invariant probability concentrates around a stable point of 
the ODE u' = b{u) (the point 3), and the more coefficient of diffusion increases, the more the invariant measure 
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Figure 3. Approximation of the stationary density for different values of c. 



is spread out. For c — 0.75 we show that the invariant measure is infinite at 0, and for c = 1 the invariant 
measures seems to be the Dirac mass at 0. 



4. Numerical example in dimension 2 

To conclude this paper we give a numerical example in dimension 2. We represent the empirical measures of 
the Euler scheme in a degenerate situation where there are two invariant measures. The first one is the Dirac 
mass at A = (0, 0) and the second one is a probability measure on R^\ {(0, 0)}. 

We consider the deterministic Van der Pol equation defined by 



y, 

(1 — x^)y — X. 



(29) 



This non-linear system of R^ has a stable point (0, 0) and an attractive limit cycle. If we add sufficient strong 
noise in whole space but not in (0, 0), the point (0, 0) becomes an attractive point for the stochastic system and 
the invariant measure of this system is ^(o,o)- 
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Write b{x,y) = ( /, 2\ I ^'^'^ '^ix,y) = ( n )■ We consider the following perturbed Van dcr Pol 

equation 

dut = b{ut)dt + (7{ut)dBt 
with ut = {Xt,Yt) e R^. We discretize the solution (X^)^^ using a decreasing step Euler scheme (X„)^^ 
and the step sequence (7„) is defined by 7„ = 0.5n~^'^. To guarantee stability of this scheme, we replace 

the function b by the function b{x,y) = I /-, ^ a a\ ) ■ "^^^ scheme is thus defined by Xq — (1, 1) and 

for every n ^ 

Xn+l = Xn+ 7„6(X„) + y/^a{Xn)Un+l, 

where [/„ is a normalized Gaussian of R^ . 

An approximation of the density of v is done using an histogram of j/^ for n — 10^ and ?7„ = 1. The histogram 
is built using a step h = 0.2. The results are presented in Figure 01 

We note that for a small noise (c = 0.5) the point (0, 0) seems not charged by the invariant probability. And 
for higher values of c the invariant probability concentrates in a neighborhood of (0, 0). 
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Figure 4. Approximation of the stationary density of the perturbed Van der Pol equation 
for different values of c. 



